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METHOD OF CALIBRATING A SPECTROSCOPIC DEVICE 

FIELD OF THE INVENTION 

ITie present invention relates to the fields of spectroscopy, 
5 spectrophotometry, and chemometrics. In particular, the present 

invention relates to a method of calibrating a spectroscopic device for 
predicting analyte levels in a non-invasive manner. The method of the 
present invention is particularly suitable for blood glucose prediction 
based on near-infrared spectrophotometry measurements. 

10 BACKGROUND OF THE INVENTION 

Biotechnological analysis and examination are often based on the 
measurement of various chemical analytes in the composition of a 
biological matrix such as blood, interstitial fluid, or living tissue. Such 
measurements may be used to evaluate a patient's state of health and to 

15 determine what, if any, treatment is necessary. For example, the 

frequent monitoring of blood glucose levels in diabetic persons witti 
glucose meters is often necessary to allow such persons to manage the 
diabetes mellitus disease, by taking insulin injections or oral dmgs to 
lower blood glucose when reqtiired. Intensive treatment based on 

20 frequent blood glucose measiirements can significantly reduce the 

incidence of blindness, kidney loss, and other complications associated 
with diabetes. 

Most home-based glucose measxirement systems require the 
patimt to invasively collect a blood sample, by pricking his or her finger, 

25 placing the sample on an appropriate test strip, and then testing the 

sample in an optical glucose meter. For millions of diabetics around the 
world, the use of lancets or otiier sharp instruments to draw blood for 
monitoring their insulin levels is a painful process, and one that often 
builds up calluses on fingers, making the collection of blood even more 

30 difficult This invasive procedure may be especially difficult to perform 

on children and flierefore particulcurly trying on parents. Furthermore, 
the test strips required for each blood sample are generally not reusable, 
and when multiple measurements are taken each day, aimoimt to 
significant costs from tite patient* s point of view. Thus, despite the fact 

35 that a large number of diabetics should take several measurements 
diroug^out each day (for some individuals, physicians recommend 
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testing glucose levels from 4 to 7 times daily), due to tiie pain, cost, and 
inconvenience involved, many diabetics do not monitor their glucose 
levels frequently enough. 

Spectroscopy is based on the analysis of how incident radiation 
interacts with the vibrational and rotational states of molecules which 
are of ahaljrtical interest. Spectroscopic measurement techniques have 
gained increased popularity because of the ability to provide fast and 
non-invasive measurements of concentrations of different chemicals or 
anal)rtes. For the reasons indicated above, this is particularly desirable 
for home based glucose meters. Spectrophotometry is a type of 
spectroscopy commonly used to quantitatively measure concentrations 
of analytes based on spectral energy distribution in the absorption 
spectrum of a sample solution or medium. In spectrophotometry, the 
energy distribution is typically analyzed within a range of die visible, 
ultraviolet, infrared, or near-infrared spectra. For example, near-infrared 
radiation (NIR) is electromagnetic radiation having a wavelength of 
between about 0.75 and 2.5 micrometers (i.e. from 150 to 400 THz). 
Near-infrared spectrophotometry generally uses instruments with 
quartz prisms in monochromators and with lead sulfide photoconductor 
cells as detectors to observe absorption bands, and NIR 
spectrophotometry is increasingly being used to measure in vivo 
analytes such as glucose, total hemoglobin, albiimin, cholesterol, and 
ethanol. 

Non-invasive, spectrophotometric meast^ement of glucose in 
human beings is performed by focusing an incident radiation source (or 
sources) on a specific part of Ihe body and detecting the spectral 
distribution of the radiation transmitted therethrough. The absorbance 
of light from the incident radiation is due to the chemical components 
within that body part including water, fat, protein, hemoglobin, melanin, 
and glucose. One difficulty with glucose measurement spectral analysis, 
is the spectral overlap between glucose and other chemicals found in 
blood, often in much greater quantities than glucose. In addition, the 
thickness, color, and structure of the skin, bones, and blood through 
which the incident radiation passes will affect the transmission (or 
reflection/absorption) tixereof. Furthermore, tite concentration of 
analytes can vary with changes in activity level, diet, hormone 
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fluctuations, and other factors. Glucose concentration measurements are 
also particularly susceptible to variations in physical and chemical 
conditions including temperature, presstire, humidity, and skin 
hydration. As a result, to perform a reliable non-invasive glucose 
5 prediction, NIR spectral measurements should be performed through a 
vascular equilibrated region of the body, and a NIR glucose 
spectrophotometer must be carefully designed so that tiie quality of raw 
spectral information from an NIR glucose meter is high. See generally 
Waynant and Chenault, "Overviev/ of Non-Invasive Optical Glucose 

10 Monitoring Techniques", IEEE LEOS Newsletter, vol. 12, no. 2 (April 

1998); and Burmeister and Amold, "Spectroscopic Considerations for 
Noninvasive Blood Glucose Measurements witii Near Infrared 
Spectroscopy", IEEE LEOS Newsletter, vol. 12, no. 2 (April 1998). 

Near-infrared glucose measurements are generally suitable for 

15 tissue depths ranging from about 1 mm to 10 cm, and are often 

performed through a patient's finger tip, although other areas of the 
body (for example the web between two fingers, an ear lobe, or tihe 
upper lip) can also be used. The sample thickness is an important 
experimental parameter because a greater thickness increases the 

20 amount of absorption and thereby lowers the minimum limits for 

detection, whereas because less incident radiation successfully traverses 
through a thicker sample (i.e. wifliout being absorbed) effectively 
increasing the spectral noise: see Bxirmeister and Arnold, supra. 

In addition to the difficulties of obtaining accurate spectral 

25 measurements with NIR spectrophotometry, a more significant 

difficulty associated with such measurements has been the need to 
calibrate such an instrument so tiiat it may be used by various different 
individuals, whose analyte concentrations and variability, for example, 
may differ substantially. Calibration of spectrophotometers (and of 

30 analytical instruments in general) is necessary to ensiure the accuracy of 

measurements performed by such devices. 

Two approaches to calibrate and predict analyte (e.g. glucose) 
levels have been \ised in the prior art In the first approach, a thorough 
calibration process is performed separately on each individual witii 

35 whom the instrument is to be used. This individual calibration process 
requires taking a number of blood samples over a period of time from 
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the individual, and obtaining reference glucose concentration 
measurements from these samples. A number of corresponding non- 
invasive spectroscopic measurements are taken concurrently, and 
calibration regression anal)^is is then performed to correlate, e.g. 
5 through linear regression analysis, the non-invasive measurements with 

the reference measurements. In this maimer, the spectroscopic device or 
spectrophotometer is custom calibrated to the specific individual. A 
major disadvantage of this approach, however, is that the calibration 
model is limited to data from that particular individual which may have 

10 an insufficient amoimt of information regarding variation of spectra 

overlapping (interfering witii) those of the analyte of interest. It is 
important to be able to characterize such spectra to enable hig^y 
accurate analyte predictions. To alleviate this problem, the custom 
calibration process may have to be carried out frequently, even on a 

1 5 daily basis. 

An alternative calibration approach requires tiiat a imiversal or 
general calibration algorithm be developed. This approach is based on 
the notion that by making a large nxmiber of calibration measurements 
(as in the first approach) a single calibration model that sufficiently 

20 accoxmts for all individual variability of all overlapping or interfering 

anal3rtes can be calculated. Chemometrics, the application of 
mathematical, statistical and formal logic methods to chemistry, is 
generally used to process and compute the spectral intaisity data and to 
produce a calibration model: see Small and Arnold, 'T)ata Handling 

2 5 Issues for Near-Infrared Glucose Measurements", IEEE LEOS Newsletter, 
vol. 12, no. 2 (April 1998); and Shaffer, Small, and Arnold, "Genetic 
Algoritiun-Based Protocol for Coupling Digital Filtering and Partial 
Least-Squares Regression: Application to the Near-Infrared Analysis of 
Glucose in Biological Matrices," Anal. Chem. 68, 2663-2675 (1996). 

30 However, for measurements of glucose, cind other analytes which vary 

distinctly and differently from one individual sample to another, a single 
universal calibration algorithm is often ineffective and can result in 
significant, and sometimes dangeroxis, erroneoiis predictions of analyte 
concentration levels. 

35 Thus, while NIR meastirements provide a non-invsisive, fast, 

painless, and convenient technique to monitor glucose levels, correlation 
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and dinical interpretation of spedial measurements to obtain the true 
glucose levels is crucial for proper therapy and disease management. 
Proper calibration of instruments for different patient populations 
(which win vary in ethnicity, age, weight, and so on) is crudal in 
5 obtaining accurate glucose prediction models. Furthermore, careful 

validation and testing of the non-invasive results and the glucose 
prediction equation is needed to determine if the glucose correlation is 
consistent in all clinically important conditions and for all, or at least 
most, t3rpes of patients. 

10 SUMMARY OF THE INVENTION 

In a first aspect the present invention provides a method of 
calibrating a spectroscopic device for providing a non-invasive 
measurement of an analyte level in a sample, comprising: 

(a) providing.a plurality of calibration algorittuns; 

15 (b) taking a set of non-invasive measurements on said sample 

with said spectroscopic device; 

(c) calculating a predicted set of analj^e levels for each of the 
calibration algorithms in response to the set of non-invasive 
measurements, each of the predicted sets of analyte levels being 

20 characterized by a variability range, a slope, an (a square of the 

correlation between said set of non-invasive measurements and said 
predicted set of analyte levels), and a standard error of prediction; and 

(d) selecting an appropriate calibration algorithm by using a 
suitability score bcised on the variability range, the slope, the and the 

25 standard error of prediction for each of the predicted sets of analyte 

levels. 

Preferably, step (d) comprises a method according to claim 1 
wherein, step (d) comprises 

(i) selecting the predicted sets of analyte levels in which the 
30 standard error of prediction is less than an upper error limit, the 

variability range is greater than a lower range limit, and the slope is 
between a lower slope limit and an upper slope limit said lower and 
upper slope limits defining an acceptable slope range; 

(ii) for each of the predicted sets selected in step (i), calculating 
35 a suitability score in response to the slope, the R^ and the standard error 

of prediction for that predicted set, and selecting the calibration 
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algorithm corresponding to tiie predicted set having the optimal 
(highest) suitability score as the appropriate calibration algorithm; 

(iii) if no predicted sets are selected in step (i), selecting the 
predicted sets of analyte levels in which the variability range is lower 

5 than the lower range limit and in which the standard error of prediction 

is less than the upper error limit, 

(iv) from each of the predicted sets selected in step (iii), selecting 
the calibration algorithm corresponding to the predicted set having the 
lowest standard error of prediction as the appropriate calibration 

1 0 algorithm; and 

(v) if no predicted sets are selected in step (i), or step (iii), 
determining that no calibration algorithm is appropriate. 

Also preferably, the acceptable slope range may be subdivided 
into a plurality of ubranges corresponding to a plurality of levels, 
1 5 comprising a first level and one or more subsequent levels, and if no 

predicted sets are selected in the first level, then repeating step (i) at each 
subsequent level xmtfl a predicted set is selected or there are no more 
. subsequent levels. 

Preferably, tiie plxirality of calibration algorithms provided in step 
20 (a) are first generated by the steps of: 

(i) compiling non-invasive and corresponding reference data 
sets of measurements of analyte levels for each of a nimiber of samples; 

(ii) rqecting^data sets tiiat are not suitable for calibration; 

(iii) combining data sets that are suitable for calibration into a 
25 plurality of groups depending on whether correlations of the combined 

data sets meet predetermined criteria; and 

(iv) generating a calibration algorithm for each of the groups 
of data sets. 

Steps (a), (c), and (d) of the method of the present invention may 
30 be performed on a computer associated witii the spectroscopic device. 

Furthermore, the sample may be an individual patient, the spectroscopic 
device may be a near-infreoed spectrophotometer, and the analyte may 
be selected from the group consisting of glucose, hemoglobin, albimun, 
cholesterol, and ethanoL 
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Hie objects and advantages of tiie present invention will be better 
understood and more readily apparent with reference to ttie remainder 
of the description in conjimction with the accompanying drawings. 
BRIEF DESCRIPTION OF THE DRAWINGS 

In the drawings which illustrate^ by way of example, a preferred 
embodiment of the invention: 

Figure 1 is a correlation scatterplot of glucose predictions for a 
number of different patients based on a universal, single calibration 
model; 

Figure 2 is a correlation scatterplot of glucose predictions for a 
number of patients based on the multiple algorithm calibration model 
according to the present invention; 

Figure 3 is another correlation scatterplot of glucose predictions 
for a number of patients based on the multiple algorithm calibration 
model according to the present invention; 

Figure 4 is a correlation scatterplot of glucose predictions for a 
typical patient based on a universal, single calibration model; and 

Figure 5 is a correlation scatterplot of glucose predictions for the 
same typical patient based on the multiple algorithm calibration model 
according to the present invention. 
DETAILED DESCRIPTION OF THE INVET^ON 

As mentioned, the present invention relates to a method of 
calibrating a spectroscopic device for predicting analyte levels in a non- 
invasive manner. The present invention can be used with a typical NIR 
spectrophotometer system having a ligjit source which is projected 
through the item to be examined, a sample interface mechanism, a 
spectrometer to separate the light into its component wavelengths, a 
detector, amplification electronics and a computer. By measuring the loss 
(absorption), between the source and the detector and applying 
appropriate chemometric (mathematical) techniques, it is possible to 
non-invasively determine the cheniicals being examined since different 
chemicals absorb different amounts of light at different wavelengths. 
Sudi a spectrophotometric device and method are described in detail in 
United States Patent No. 5,361,758, the contents of which are 
incorporated herein by virtue of this reference. While the present 
description relates primarily to glucose measurement, one of tiie major 
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fields of application for NIR measxirement at present, it will be 
understood that the principles of ttie present invention eqtdvalently 
apply to other analytes non-invasively measured using various 
spectroscopic techniques. 

To utilize the NIR spectrum for glucose measturement, it is 
necessary to use a spectrometer which has wide dynamic range, a high 
signal to noise ratio, and e>diibits low scattering losses. The output from 
the spectrometer is used to generate q>ectra v^tti high precision both in 
absorbance and wavelength. Significant glucose absorption bands are 
centered about wavelengths of 1.67, 2.13, 2.27, and 2.33 micrometers (as 
discussed in Small and Arnold, "Data handling Issues for Near-Infrared 
Glucose Measurements", supra). In addition to these wavelengths there 
are significant glucose absorption bands at 0.97 and 1.12 micrometers 
which wavelengths allow transmission measurement through a greater 
distance in the tissue. To be able to use NIR to measure a particular 
compoimd/ ancQyte, chemometric mathematical analysis is applied to the 
measured spectrum. The mathematical analysis techniques are carried 
out by a computer equipped with advanced software capable of 
interpreting the resulting complex spectra. 

To universally calibrate an NIR spectrophotometer, the 
absorbance associated with the compoiind or anal3rte of interest must 
first be measured on a relatively large number of samples. These NIR 
measurements are then compared to measiurements made in a more 
traditional and more accurate invasive manner. From these comparisons 
an algorithm is developed that characterizes the analyte to be measured. 
The methods used to generate calibration models for NIR 
spectrophotometers are sophisticated- An important criterion is that 
calibration samples be uncorrelated witii other chemical spedes 
concentrations. Satisfaction of these criteria is burdensome because of 
the large nimiber of chemical species typically present in a biological 
matrix or growth medium. 

The present invention provides a multi-algorithm calibration 
method for use with a spectroscopic device capable of predicting analyte 
levels in a biological matrix, such as blood glucose levels, for a variety of 
different samples with significantly greater accuraqr than prior art 
tiniversal single calibration models. The calibration meftod of the 



wo 01/16579 



PCT/CAOOA)1005 



- 9 - 

invention has two key parts: the generation of a plurality of different 
algorithms and the selection of the algorithm which is most appropriate 
for generating an accurate prediction of analyte concentration. 

To generate suitable calibration algorithms, the non-invasive and 
reference measiarement data sets for a number of individuals are fiist 
compiled. Once this is achieved, correlation techniques are carried out to 
relate non-invasive spectral measurements to reference glucose levels. 
At fliis stage the suitability of individual data samples for calibration 
purposes must be considered. If some of tiie non-invasive q>ectral 
measurements in a data set resvdt in a saturation of the measurement 
equipment, e.g. at the output of an analog-to-digital converter, then the 
data sample should be disqualified for use in a calibration algorithm. 
Sinularly, if the range of variability of the non-invasive spectral 
measurements in a data set is too small, the data sample should also be 
disqualified from calibration. Furthermore, if the correlation of glucose 
in a particular data set is imdesirably correlated to another variable, such 
as the concentration of another chemical, the sample should again be 
disqualified. As a result, a large nimibers of sample sets will not be 
suitable for calibration algorithm development. 

The grouping or combining of calibration-suitable data sets (each 
specific to a particular individual) into different algorithms occurs on a 
trial and error basis, by repeatedly combining these data sets and 
performing a new correlation on the combined sets of data. If the 
correlations of the combined sets meet certain criteria, then the 
combination is acceptable, if not, a different combination is attempted. 
Once a combination is deemed acceptable, the addition of a further data 
set is attempted and its acceptability is determined in turn. In this 
manner, a combination of data sets for an algorithm grows from 2 data 
sets, to 3 data sets, to 4 data sets, and so on until a sufficiently large 
number of data sets is obtained from which a reliable calibration 
algorithm can be developed. 

In combining data sets for a single algorithm, two criteria are key. 
First, any correlation of glucose in the combined data set group with 
other parameters must be miiumized. Such imdesirable correlations may 
be exhibited with different chemical species concentrations or other time 
varying quantities. A suitable quantitative test for ihis criteria is to 
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ensiire that all undesirable correlations of the combined data sets, as 
indicated by the square of the correlation coefficient R^^ (for the 

undesirable correlations), are less than a certain limit, such as R^^ < 0.10. 
In known manner, the correlation coefficient, R^, is a measure of the 

5 tendency of two variables to vary together. If the addition of a data set 

to another data set or to a suitable sub-combination of data sets does not 
satisfy this criteria, tiie new combination will not be acceptable. 

In addition, the combination of data sets can be optimized by 
ensuring that the data sets exhibit a high degree of correlation with one 

1 0 another, so a second criteria in combining data sets should require that 
the data sets in a particular group be sufficiently correlated with one 
another. This can be achieved, for example, by ensuring that the 
Standard Error of Prediction (SEP) is minimal and the square of the 
correlation coefficient, R^, together with tiie slope relating reference 

1 5 glucose to predicted glucose, are greater than about 0.8 for an existing 

sub-combination of data sets and another data set being potentially 
combined with that sulxrombination. 

In an alternate embodiment of the present invention, iiostead of 
successively adding on to, or building up, groups of data sets for a 

20 calibration algorithm, a large group of calibration suitable data sets can 

be initially taken up at random. In this embodiment, data sets are 
removed from the group, and the SEP and correlation criteria above are 
reevaluated. If tiie SEP, slope and correlation criteria improve after the 
removal of the data set, the data set is permanently removed feom the 

25 group, whereas if the criteria do not improve, the data set is reinstated 

back into ttie group. Using this top-down approach, the removal of data 
sets can cease when the criteria meet accq>table limits. 

The SEP, slope/ and correlation calculatior\s for determination of 
the calibration algorithm groupings (and the associated calibration 

30 zdgorithms) become increasingly complex with increasing amounts of 

data and are preferably carried out using partial least-squares (PIS) 
regression analysis. The PLS analysis techniques are well known to those 
skilled in the art and provide good linear approximation by removing 
considerable redundant infonnatioru see generally Small and Arnold, 

35 'Data handling Issues for Near-Infrared Glucose Measurements", supra. 
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Other types of anal)^is, such as partial component analysis (PCA) or 
artificial neural networks (ANN), can also be used. 

This calibration algorithm generation process is used to generate 
at least two groups of data sets (corresponding to groups of individiials), 
5 each providing a calibration model or algoritiim for predicting the level 

of anatyte present in a biological matrix based on a non-invasive 
spectroscopic measurement Preferably, two initial groups of data set are 
developed, and, from these two groups, additional groups can be 
generated. These additional groups are generally smaller in size tiian the 

1 0 two initial groups and may meet more stringent correlation criteria. 

These additional groups may be constructed with tiie top-down 
approach described above, wherein data sets are selectively removed 
from an initial group to determine if the correlation criteria improve or 
deteriorate after the removal. 

15 Given the availability of multiple calibration algorithms, it is 

necessary to select the most suitable of these algorithms for predicting 
the levels of analyte in any given patient or sample. Once this is done the 
spectroscopic device is calibrated with the selected algorithm. The 
algorithm selection process is generally independent of the process used 

20 to generate the multiple algorithms, altiiough it will dearly depend on 

flie number of different algorithms that are available. 

For any given patient or sample on which the non-invasive 
spectroscopic measurement is to be made, it is necessary first to compile 
a set of highly accurate reference measurements and corresponding 

25 non-invasive spectral measiu'ements at selected intervals over a period 
of time. In a preferred embodiment, two successive non-invasive 
spectral measurements are successively taken each measurement 
interval, and their mean is then taken as a single non-invasive 
measurement for that interval. 

30 Before algorithm selection begins, two variables may be set to 

determine the stringency of acceptance requirements for patients (i.e. 

whether at least one of the available calibration algorithms wiU be foxmd 

suitable for a patient). These are (i) a lower range limit for the range of 

glucose or analyte variability in the data set (if the data does not vary 

35 significantiy it may not provide sufficient information to be useful), and 

(ii) an upper error limit (SEPj,^) for the standard error of prediction 
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(SEP), ie. the square root of the mean of the squaied deviations of the 
reference analyte concentrations from the anal3rte concentration values 
predicted by a calibration algorithm. 

From ttie non-invasive spectral measurements, a predicted set of 
analyte concentration levels are calculated according to each calibration 
algorithm. When compared to the reference set of concentration levels, 
each predicted set will be calculated from a linear regression or 'Tjest fit" 
line characterized by: a range of glucose level variability; a slope (which 
ideally equals 1, i.e. the predicted values correspond directly to the 
reference values); a correlation coefficient R, defining the correlation of 
tiie predicted value set with the reference value set; and a standard error 
of prediction (SEP), as defined above. Note that as the correlation 
coefficient R increases the SEP will decrease, and vice versa, exc^t in 
cases in which there is an offeet or bias in the predicted analyte. 

By way of example, and not to restrict the scope of the present 
invention, the algorithm selection process can include 4 levels or steps. In 
a first level, if the predicted data set for an algorithm meets the following 
criteria 

i) SEP < upper error limit (SEPji^) 

ii) Analyte range > lower range limit 
and iii) 0.5 < slope (of regression line) < 1.05, 

a suitability score is calculated according to the following equatioru 

score = (slope) (R^) / (SEP) 

If more than one algorithm meets the first level criteria, die algorithm 
wiHx the highest smtability score is selected and the selection process is 
complete. If no algoritiun qualifies in terms of the first level criteria, die 
selection process advances to a second level. 

The second level is identical to the first level, except that the slope 
criteria becomes 
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For algorithms that meet this revised slope criterion (and the SEP and 
range criteria), a suitability score is calculated as in the first level. Again, 
if more than one algorithm qualifies at the second level, the algorithm 
with the highest score is selected, and the selection process is complete. 
5 At a third level, if none of the algorithms qualify at either the first 

level or the second level and the range of anal)rte measxirements is 
greater than the lower range limit, the sample or patient is excluded 
from measurement, since under these circumstances no calibration 
algorithm can be relied on to perform accurately. 

10 Lastiy, at a f oujih level, if none of the algorithms qualify at either 

the first level or the second level and the range of analyte measurements 
is lower than the lower range limit, the algorithm with the lowest SEP is 
selectedif that SEP is less than the upper error limit If the SEP is greater 
than ttie upper error limit, the patient is again excluded from 

1 5 measurement. Thus, in the fourth level, the SEP becomes the sole factor 

in deciding whether an algorithm should be accepted. 

Note tfiat the slope criteria limits in the first and second levels 
described above are preferred limits, and these levels can also be varied 
somewhat (similar to tfie lower range limit and upper error limit) 

20 depending on the application and acceptance requirements. Also, the 

number of levels or steps used in the algoritiim selection process can 
vary, from a single level covering the entire preferred range, to two or 
more levels subdividing the preferred range into two or more range 
limits or range criterion. Regardless of the number of levels, however, 

25 the third and fourtii levels as described above will correspond to the last 

two levels, in which the patient is either excluded, or a decision to accept 
an algorithm is made on the basis of SEP alone. It will be readily 
appreciated by those skilled in the art that ttie number of levels or steps 
chosen is a matter of preference only. 

30 To illustrate the principals of the present invention, the multiple 

algorithm calibration method was used to calibrate a NIR 
sj>ectrophotometer for measuring glucose levels in a number of 
different patients. AU reference data measurements were taken with a 
glucose anal3^zer such as the industry standard YSI™ (Yellow Springs 

35 Instruments) Glucose Analyzer which provides precise glucose 

measurements on (invasive) blood samples. The inventor foimd tfiat 
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only approximately 6% of single NIR patient data sets qualified for 
calibration purposes, given the criteria described above. These 
calibration-qualified sets of data were grouped into two initial calibration 
algorithm sets, each initial group was built up to a group with 33 sets of 
data within it Eleven additional groups were spun off from these two 
initial groups and each of these groups had at least 8 sets of data within 
it. Thus, tiiirteen different calibration algorithms were available for 
predicting glucose levels. 

The algorithm selection was based on a total of 24 reference 
measurements (i.e. YSI™ measured) and 24 means of two NIR finger 
measurements, taken from a patient over a three day period. The above 
described algorithm selection proc^ was run for a nxunber of ditferent 
patients (not involved in the calibration process). 

Figtire 1 shows a correlation scatterplot of glucose predictions for 
a number of different patients based on a imiversal calibration model 
(QUAD 37). The scatterplot in Figure 1 is superimposed on a Clarke 
error grid, as described in Clarke and Cox, "Error Grid Analysis", 
Diabetes Care, 10:622-628 (1987). The Clarke error grid breaks the 
correlation space into five regions (A-E) that assess measurement 
accuracy on the basis of validity of the corresponding clinical dedsiort 
Correlation points falling within the "A" region correspond to ttie correct 
clinical decision being made based on the similarity between the actual 
and predicted glucose levels. (Le. the predicted values deviate by no 
more than 20% from tiie reference values). In the "B" region, the 
predicted values deviate by moie ttian 20% from the reference values, 
but treatment decisions made based on the predicted levels of glucose 
would not jeopardize or adversely affect the patient In regions "C", "D", 
and "E", the predictions significantly deviate from the reference values, 
and treatment decisions based on these predictions may well be harmful 
to a patient. The Clarke error grid is often used to evaluate the diiucal 
consequences of home-use blood glucose monitor errors, in contrast to 
more conventional analjrtical methods which may be more suitable for 
laboratory reference devices. 

Referring to Figure 1, it can be seen that the upper error limit 
S^max imiversal calibration algorithm is 3.07 mmol/L and tt\at 

significant amounts of the data are outside acceptable limits or regions. 
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In contrast. Figure 2 shows the correlation scatterplot of glucose 
predictions for 25 patients based on the multiple algorithm calibration 
model according to the present invention, witfi Stir'^ax ^q^^ to 2.76 

mmol/L. Similarly, Figure 3 shows a correlation scatterplot of glucose 

5 predictions for 22 patients based on the multiple algorithm calibration 

model according to the present invention, witii SEPj^^ equal to 22B 

mmol/L. Both the predictions in Figure 2 and particularly Figure 3 
exhibit much better accuracy than those in Figure 1 (calibrated according 
to a imiversal algoritiim). 

1 0 Similarly, Figure 4 shows a set of correlation data for a typical 

patient, calibrated according to ti\e same universal algorithm, whereas 
Figure 5 shows the correlation data for tiie identical patient calibrated 
according to the mtdtiple algorithm method of the present invention. 
Once more, the improved accuracy of the glucose predictions based on 

1 5 the calibration method of the present invention is evident. 

It will be dear to those skilled in the art that the principles of the 
present invention are applicable to a broad range of spectroscopic 
applications. For example, NIR radiation can be used to penetrate such 
items as hximan tissue, vials of blood, or containers of milk - all normally 

20 considered opaque. Mathematical analysis of the resulting absorption 

spectrum determines the composition of the substance penetrated by 
the light. Thus, in a blood testing lab, samples can be examined without 
contact by passing the light through a plastic vial containing ttie liquid. 
Similarly, in a dairy application, the NIR light is passed through a bag of 

25 milk and such parameters as butterfat, solid particles and lactate are 

measured. The multiple algorithm method of the present invention is 
easily extendible to these types of non-invasive measurements as well. 
Furthermore, the method of flie present invention can also be applied to 
the non invasive measurement of a mmiber of otiier blood analytes sud\ 

30 as cholesterol/ hemoglobin, HblAc, Fructosamine, and 1.5 AG. 

While preferred embodiments of the present invention have been 
described, the embodiments disclosed are illustrative and not restrictive, 
and the invention is intended to be defined by the appended claims. 



